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' Using the complex energy method, the problem of nucleon-deuteron scattering is solved 

I with a simple three-body force having a separable form. Our results are compared with the 

results of modern direct two- variable calculations and a good agreement is found. This forms 
a firm base for other applications of the complex energy method. 

§1. Introduction 

^ ■ Faddeev calculations of nucleon-deuteron (Nd) scattering with a three-nucleon 

■ force (3NF) have been performed for many years.l^ Modern nucleon-nucleon poten- 

tial stem from the chiral effective field theory.^^ This theory gives us also various 
3NFs but for many years numerous calculations were performed with the Tucson- 
>'. Melbourne 3NFP There are two alternative formulations leading to the scattering 

\ amplitude with 3NF in the momentum space.® One of them was introduced by 

00 ■ Glockle and Brandenburg.'^ In this scheme the Alt-Grassberger-Sandhas (AGS) 

. equation for the elastic scattering transition operator U is written a^ 

(N 

O: U = PGo^ + {I + P)t^ + PtiGoU + {1 + P)tiGotiGoU, (M) 

(N- 

where P, ti, Gq and t4 are the permutation operator (P = P12P23 + -P13-P23)) the 
two-body t-matrix, the free Green's function and the t-matrix generated by a 3NF, 
j/J ; respectively. 
'V^ ■ The other approach has been used by the Bochum-Krakow grou (BK). Here 

5^ . the equation for U has the following form: 

U = PGo^ + (1 + P)vi^\l + P) + PtiGoU + (1 + P)vi^\l + P)GotiGoU, (1-2) 
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where V4 is a part of 3NF which is symmetric under the exchange of nucleons 2 

and 3 so the full 3NF, V4, is given as V4 = V"4^^(l + P). In practice, first the auxiliary 
equation for the three-body operator T is solved, 

T\^) = hP\cl)) + (1 + hGo)V^\^) + hGoPT\^) + (1 + tiGo)V^GoT\^) , (1-3) 

with the incoming momentum space state \(j)) composed of the relative nucleon- 
deuteron motion and the deuteron wave function. Using this T, the elastic transition 
operator U is given as 

U\^) = PGo V) + PT\^) + V^\<P) + ViGoT\cl)). (1-4) 

The consequence of employing Eq. (|l-3p is that the t-matrix of the 3NF, ^4, does not 
appear in this formulation. The formalism of Eq. (jl-ip has already been used with 
a simple separable 3NF.'^ In that calculation, however, the t-matrix of the 3NF was 
approximated by taking t4 = V4. It is the aim of the present investigation to show 
that both schemes, (jl-ip and (|l-3p - ()l-4p . yield the same numerical results under a 
consistent separable approximation for the underlying dynamics. 

In Section 2 we show details of our formalism which leads to an analytical t- 
matrix of 3NF and the modified Amado-Lovelace (AL) equation. In Section 3, using 
the Yamaguchi separable potential and a separable 3NF, we will demonstrate results 
for the Nd elastic scattering cross section and the neutron-deuteron scattering length. 
These results will be compared with the results of standard calculations using the 
general two-variable formulation.^ In Section 4 we give a short summary. 

§2. Modified Amado-Lovelace equation 

We first rewrite Eq. lfTT]) aP 

U = hGoPh + tiGo(l + P)UGQti + tiGoPU + tiGo(l + P)UGqU , (2-1) 

where U is related to U via U = tiGoUGoti and the t-matrix of the 3NF t4 is the 
solution of the Lippmann-Schwinger equation: = V4 + ^4^0*4. Next we introduce 
the so-called separable approximation for the two-body and three-body interactions. 
We choose first the Yamaguchi for mf actor 9c{p), for the pure S-wave two-body 
potential: 

Vc{p,p') = -Xc9c{p)gc{p') = -Ac ^2|^2 (2-2) 

where p' (p) is the magnitude of the initial (final) relative momentum within the 
two-nucleon subsystem and Ac is the coupling coefficient. We choose the same value 
for the /3 parameter (1.4488 fm~^) both for the ^Sq (c = (j)) and ^Si (c = d) states. 
The coupling coefficient A^ for the ^Si state is taken as A^ = 87r/3(a + /3)^/m to 
get the experimental deuteron binding energy with a =0.2316 fm~^ , and the other 
coupling coefficient X^p = 87r/3(ao — (5)"^ /m for the ^Sq state allows us to reproduce 
the nucleon-nucleon phase shift parameters with oq =-0.040 fm~^. The Lippmann- 
Schwinger equation, ti = v + vG^ti , is explicitly given as: 



A Benchmark Calculation for the Nd Scattering with a Model Three-Body Force 3 



t,ip,P;E) = v{p,p) + J^ + 2,2 ■ (2-3) 

where m is the nucleoli mass (1/m =41.47 MeV fm^). It is well knowiP that the 
t-matrix ti of a separable potential can be calculated analytically and has also a 
separable form {ti{p,p'; E) = g{jp)T{E)g{p')). 

In the case of the 3NF we assume the following S-wave phenomenological sepa- 
rable potential 

V,{p,q,p'ci') = XMp,<iMp',<i') = ^4 + 3^2 + ^2)2 (^,2 + 3^. + ^2)2' (2-4) 

where q is the Jacobi momentum for the spectator particle and A4 is the 3NF coupling 
coefficient. The choice of the A parameter will be discussed in Section 3. Like the 
two-body t-matrix, also the 3NF t-matrix t4 can be obtained analytically, again in a 
separable form: t4{p,q,p' ,q'; E) = h{p,q)Ti{E)h{p' ,q'). The derivation is presented 
in the Appendix. 

Using these separable formulas for ti and t^ leads to a modified AL equation: 

Xc,c'{q,qQ]E) = Zc,c'{q,qo]E) 

poo o //2 n"'^tin" 

+ E I z,,c"{q,q";E)T,„{E - ^)x,„,Aq",qo;E)^^. (2-5) 

The amplitudes X^^c' and Z^^c' are written as 

(p, g, c\U\p', q', d) = gc{p)Tc{E - ^)Xc,c'{q, q'; E)tc'{E - -^)gc'{p), (2-6) 

{p, q, c\tiGoPti + tiGo(l + P)tiGoti \p', q', c') 

= 9c{p)tc{E - —)Zc,c'{q,q';E)Tc'{E - -^)9c'{p') 

g,{p)r,iE - ^){zj^%,q'-,E) + Zi%,q'-,E)}TAE - ^)gAp') (2-7) 



with 



and 



J_i E — [q^ -\- q''' -\- qq'x)/m + le 



7(3)/ /. _ r 9c{p)h{p,q) p dp 

Z^^,{q,q^E)- ^_ +3^2/4)/^ + ,,^ ^^i^) 

r Kp',q')gAp') P"dp', , . .... 

where Pc{x) is the Legendre polynomial. Here A^c' is the spin and isospin recoupling 
coefficient defined in the Appendix. We would like to emphasize that the quantities 
and can be calculated analytically (see the Appendix for details). If the 

/ON _____ 

Z^J term is omitted, Eq. (|2-5|) is brought back to the original AL form. 
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§3. Cross Section and Scattering Length of Nd elastic scattering 



Here we would like to demonstrate some results. The triton binding energy 
calculated with the Yamaguchi potential^ is -13.0 MeV. Comparing this number 
with the experimental result (-8.48 MeV) clearly shows overbinding. This is because 
the Yamaguchi potential is too simple and furthermore restricted to act only in 
S-waves. Corresponding numbers obtained with modern potentials show in turn 
underbinding by about 1 MeV. Since our aim is to check consistency between the 
numerical results obtained with Eq. (jl-ip and Eqs. (|l-2p - (|l-4p . we introduce a simple 
3NF as given in Eq. ()2-4p . That 3NF contains two parameters (A4 and A), which 
are chosen to reproduce the experimental triton binding energy 8.48 MeV. Figure [T] 
shows the relation between A4 and A under this condition. 





Fig. 1. The relation between the parameters 
A and A4 from Eq. (|2-4p . which is ob- 
tained by fixing the three-nucleon bind- 
ing energy to Et,^ -8.48 MeV. 



Fig. 2. The doublet scattering length as 
a function of the parameter A. Note that 
A4 and A are now correlated as shown in 
Fig. \T\ The experimental value of is 
0.65 ±0.04 fm.lS 



Before introducing 3NF our result for the doublet scattering length is = -2.17 
fm, which is very different from the experimental data, 0.65ib0.04 fm.® In the case 
of the quartet scattering length our prediction is = 6.86 fm, which is rather close 
to the data (6.35±0.02fm).9' 

Figure [2] shows results for when the 3NF is included. Though the doublet 
scattering length is sensitive to the strength of the 3NF, within the chosen range of 
the A parameter the experimental value can not be reached. However, inclusion of 
the 3NF clearly improves the situation. On the other hand, the quartet scattering 
length hardly changes with A (not shown). 

In the actual calculation we take A = 2 fm'^ and A4=16520 MeVfm^^. The 
Nd elastic scattering differential cross section is calculated at Eiab=l4:A MeV, using 
Eq. p-ip and then separately with Eqs. (|l-2p - (|l-4p . The Coulomb force is neglected 
in both our calculations. In order to achieve convergence, all 3N states with total 
angular momentum up to 31/2 for the both parities have been included. We use 
100 q Gaussian integral points and 20 x points for the angular integration. The 
maximum value of the q momentum is set to qmax= 20 fm~^. The smallest e, which 
appears in the Green's functions in the complex energy methocP is typically 0.01 
fm~^. 

In Fig. [3] the differential cross section for elastic Nd scattering is demonstrated, 
comparing calculations with and without the 3NF. The theoretical predictions ob- 
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tained with the two different schemes agree very weh with each other (up to 1 %) 
but yield only a fair description of the experimental data. 

§4. Conclusions 

The aim of the present investigation was twofold. First of all we wanted to com- 
pare two calculational schemes, which deal with Nd scattering under the inclusion of 
a 3NF. To this aim we chose a simple separable form of the nucleon-nucleon potential 
and the 3NF. In this case, in the scheme using Eq. (jl-ip . some parts of the calcula- 
tions can be done analytically, which leads to very accurate results. The comparison 
of these results with the results basing on Eqs. p-2p - (|l-4|) provides a very good 
test for the numerical performance of this second method, which treats directly any 
nucleon-nucleon and 3N forces. Secondly, since the first numerical framework uses 
the complex energy method, we could confirm that this method can be successfully 
used in the few-nucleon calculations. We plan to apply the complex energy method 
in the three-dimensional treatment of two- and three- nucleon systems.'^ 
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Fig. 3. The differential cross section of elastic Nd scattering at J5;ai,=14.1MeV. The short-dashed 
(dotted) line shows the theoretical results of the separable scheme (BK scheme) without 3NF. 
The solid (long-dashed) line demonstrates the results of the separable scheme (BK scheme) with 
the inclusion of 3NF. The cross ("x") marks are for the experimental data frorrP^ (F^)- 
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Appendix A 

Analytical expressions 

We show some analytical expressions used in the present investigation. The 3NF 
of Eq.(l23D satisfies 

U{p,q,p',q';E) = Vi{p,q,p' ,q') 
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/-oo p"2dp roo q"2^q y^^p^ ^^f^ q")U{p" , q",p', q'; E) 



JO 27r2 27r2 E - {p"'^ + ?,q"^ / 4.) / m + ie ' 

Assuming the separable form {t4{p,q,p' ,q'; E) = h{p,q)Ti{E)h{p' ,q')) we have 

Ti{E) = A4 + \iu{E)h{E) , (A-2 

where the integration /4 is performed using the so-cahed "hyper coordinates" {p^ = 
p2 + 3^2/4, tan = ^q/{2p)), 



h{E) 



00 roo 



h{p, qY p^dp q^dq 

E-{p'^ + 3g2/4)/m + ie 27r^ '2^ 



/•2/)2 /) 7/] 

sm t^cos t^dt' 



384V37r3/l2(^m + yl2)4 
X {{Em + yl2)(2£;2m2 - S^m/l^ - A"") + 6yl2£;^m2(log(-£;m) - log(yl2)))(A-3) 



where log(— ii^m) = log(ii^m) — vri in case of > 0. 
In order to calculate the Z^^J terms we use 

r P^dp h{p,q)g{p) ^ m(/3 - ik' + 2^10 

io 2^"^ E - {p'^ + 3q^/4)/m + ie 8TT{f3 - ik')A'{(3 + A')^{A' - ik')^'^ ' 

with /c' = Y^m£;-3g2/4 and yl' = y/SgVlT^. 

Finally, the recoupling coefficients A^c' for spin doublet (I = ^) and quartet 
(I = I) states in Eqs.(lM])-(l23]) are 
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